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I. INTRODUCTION 

Experimental progresses over the last few years provided exciting results in the heavy baryon sector. There have 
been plenty of heavy baryon observations yielding a large amount of experimental data on charm and bottom baryons. 
Observation of Belle, BABAR, DELPHI, CLEO, CDF, DO etc [2-i] have motivated the theoretical interests on heavy 
baryons containing c and b quarks. The mass spectros cop y of these baryons have been studied using various theoretical 
mode ls llCj| -[l8j as well as QCD sum rules method fl9[-[3ll . Their masses calculated via QCD sum rules in heavy quark 
limit [22j and in the heavy quark effective theorv [23l425j . 

Beside the mass spectrum there have been theoretical studies on the magnetic moments of these baryons utilizing 
naive quark model [32|, [33[ , quark model [HI, [3pl . bound state approach (3|| , relativistic three-quark model [37[ , hyper 
central model [38j. Chiral perturbation model |39|. soliton model skyrmion model [4l|, nonrelativistic constituent 
quark model [42j . QCD sum rules in external magnetic fiels [Hj], and light cone QCD sum rules (LCQSR) method 
0,111. 

In the present work we make use of the (LCQSR) to calculate the coupling constant g-E Q ~,Q-n- A similar work has 
been done in [46] for the coupling constant (/SqAqtt- This paper is organized as follows. In section 2, we introduce 
the interpolating field for Eq and give the details of (LCQSR) calculations for the coupling constant. In section 3 the 
numerical analysis, discussion and conclusion are presented. 



II. LIGHT CONE QCD SUM RULES FOR THE gs Q ~ Q7r COUPLING CONSTANT 

To calculate the coupling constant g-EQE.Q-K via the (LCQSR) one studies a suitably chosen correlation function of 
the form 

n = i J atxe^ (Tr(g) | T{m Q {x)m Q (0)} | 0), (1) 

where rjs Q (x) denotes the interpolating current of Sq baryon and T denotes the time ordering product. One can 
calculate this correlation function either phenomenologically, inserting a complete set of hadronic states into the 
correlator to obtain a result containing hadronic parameters, or theoretically via the operator product expansion 
(OPE) in deep Euclidean region p 2 — > 00 in terms of QCD parameters. Sum rules are obtained by matching these 
two expressions after Borel transformations and the contribution of the higher states and continuum is subtracted. 
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The calculation of the phenomenological side is similar to the calculation in |46j and the details are presented for 
completeness. To obtain the physical representation of the correlator a complete set of hadronic state having the 
quantum number of Sq baryon is inserted. Then the correlation function becomes 



n 



(Q I m Q I -q{p2)) . 

-2 — 2 (^ Q {p 2 )ir(q) | c Q (pi))- 

P2 m E 



go(pi) I m Q I Q) 

Pi ~ mi a 



(2) 



where p\ = p + q and p2 = p and ... represents the contribution of the higher states and continuum. The matrix 
elements representing the coupling of the interpolating field to the baryon state under consideration are defined as 



(0 | m Q I 2 Q (p,s)) = Xs Q m Q (p,s), 



(3) 



where Ah q denotes the coupling strength and us Q is the spinor for the Sq baryon. The coupling constant <?E Q H Q7r is 
defined by the matrix element in Eq. which is given as 

(E Q (p 2 )TT(q) | Hq(pi)) = g-E Q ~ Q ^u{p 2 )il-,u{p 1 ). 



Using the Eqs. (|3]) andQin Eq. ((2} one obtains the phenomenological side of the correlator as 



n 



gSQ5q-n-|A5 | 2 

'{p\- m| Q )(pf- to| q ) 



(4) 



(5) 



The coefficient of any one of the structures ^ ^75 or ^75 can be used. In this work, we will work with the structure 

#75- 

For the calculation of the QCD side of the correlation function which is obtained via (OPE) one needs to know the 
explicit expression of the interpolating field of Sq which is given in the following form: 



ui Cs b fQ c + (3 < C7 s h )Q C + < CQ b 7 °s c + p < C^Q b s 



(6) 



where Q represents the heavy quarks c or b, /3 is an arbitrary parameter with /3 = — 1 corresponding to the Ioffc 
current, C is the charge conjugation operator and a,b,c are the color indices. After inserting the interpolating fields 
into Eq. (JTJ) and carrying out the contractions, the following expression is obtained: 
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+Tr[ l5 CS^ a 'C l5 S bb '}S c Q c ' + Tr[ l5 CS^ a 'C l5 S b Q b ']S[ 



0). 



(7) 



Note that this is a schematical representation. The pion can be emitted from any one of the iiord quarks and hence 
both contribution should be summed. To obtain the contribution of pion emission from any one of the quarks, its 
propagator is replaced by : q a (0)q b (x) :. To proceed with the calculation the heavy and light quark propagators are 
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needed. In this work, the following propagators are used: [47 



Sq(x) = S£ ee (x)-ig s 
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where the free light and heavy quark propagators in Eq. (|8|) are given in x representation as 

ofree \ $ 

« ~~ 2ir 2 x 41 



S; 



free m Q K^mgy^X 2 ) m 2 Q j 

Q = — 7== ^^KV-x 2 ), 



4^ 2 <J 



(8) 



(9) 



where Ki are the Bessel functions. As seen from Eq. (O as well as the propagators the matrix elements of the form 
(Tt(q)\q(xi)T iq(x2)\Q) are also needed. Here I\ represents any member of the Dirac basis i.e. {1, j a , cr a p/^/2, 17570, 75}. 
In terms of the pion light cone distribution amplitudes the matrix elements (7r(g)|g(a;i)rjg(a;2)|0) are given explicitly 
asii 
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(*(p)W(x)i M ig a G e 0(vx)q(O)\O) = Pfl (p a x fi - Pp x a ) — f„rr£ j Vae l ^ +va ^ x V\\{a % ) 

px J 



+ 



i 

px 

1 



U™1 



x I Vae l(a « +va <> )px V±{a i ), (10) 

where = /„ m ^r md , Aw = > = dagda g da g 5(l - a, - a g - a ff ) and the tp^(u), A(it), B(u), ^p(u), <yS CT ("), 

7"(aj), *4_i_ (aii), ^4||(aj), V±{cti) and V||(aj) are functions of definite twist and their expressions will be given in the 
numerical analysis section. 

With these inputs, the correlation function can be calculated in terms of quark-gluon degrees of freedom. To match 
the two representation, their spectral representation is used. The contributions of the higher states and continuum 
are subtracted using quark-hadron duality. Furthermore, to eliminate the unknown polynomials in the spectral 
representation and suppress the contribution of higher states and continuum, Borel transformation is applied with 
respect to p 2 and (p + q) 2 . Finally, the sum rules is obtained from the integral: 



e^ 1 p(s)ds + e m' a T, 



(11) 



where the explicit expressions of p(s) and T are given in Appendix A. 

To obtain a prediction for gs Q s Q7r , the residue Ah q is also needed. The residue can be calculated using mass sum 
rules and is given as: 



-mV M 



el? 7 p 1 (s)ds + e « 2 T\, 



(12) 



where explicit expressions of pi(s) and Ti are given in Appendix B. From Eq. ()12|) . the mass can be obtained by 
differentiation with respect to M 2 as 



G e&s Pl {s)ds + M^ieSTT! 



Im 2 Q ^ Pi{s)ds + e~M*~Yi 



(13) 



III. NUMERICAL ANALYSIS 



In this section the numerical analysis for the coupling constant gE. Q E Q ir is presented. The required input parameters 
are given as: (uu)(l GeV) = -(0.243) 3 GeV 3 , 6(1 GeV) = 0.8(uu)(l GeV), m b = 4.7 GeV, m c = 1.23 GeV, 
ms b = 5.792 GeV, m Sc = 2.470 GeV, and m 2 (l GeV) = (0.8 ±0.2) GeV 2 [19], / w = 0.131 [Mill, m w = 0.135 GeV 
. We also need the 7r-meson wave functions for the coupling constant calculation, whose explicit forms are presented 
as H El 
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where C^{x) are the Gegenbauer polynomials, 
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The constants in the Eqs. (fl4| and (fl5|) are calculated at the renormalization scale /J, = 1 GeV 2 and are given as 
of = 0, a\ = 0.44, 773 = 0.015, 774 = 10, w 3 = -3 and w 4 = 0.2. 

Looking at the result of LCQCD sum rules calculation for the coupling constant gE Q B. Q TT one encounters three 
auxiliary parameters. These parameters are the Borel mass M 2 , the continuum threshold sq and the arbitrary 
parameter f3 and there should be no dependency of a physical quantity, such as the coupling constant for our case, 
on them. Therefore at this stage a working region of these auxiliary parameters should be determined. In order to 
determine the upper and lower bound of M we use the requirements that the continuum contribution be less than 
that of the ground state, and the highest power of 1/M 2 be less than 30°/o of the highest power of M. The former 
(latter) is used to determine upper (lower) bound of M 2 . To determine the value of continuum threshold sq we use the 
two-point correlation function from which we obtain the mass sum rules. In Fig. 1 and Fig. 2, we plot the dependence 
of our prediction on toh c to the Borel parameter M 2 , and cos6*, where j3 — cos 9, respectively. For these plots, the 
continuum threshold is chosen to be sq = 8 Gev 2 and sq — 9 Gev 2 . For these values of the continuum threshold we see 
from these figures that our predictions are in agreement with experimental results and that there is no dependence on 
the auxiliary parameters. In Figs. 3, 4, we carry out the same analysis for TO~ b and find the continuum threshold to 
be so = 36 Gev 2 and sq = 37 Gev 2 . In Fig. 4, we also observe that our prediction is stable with respect to variations 
of 9 for the region 0.5 < cos# < 0.5 which corresponds to \/3\ > 1.7. 

The results for the coupling constant calculation are presented in Figs. 5, 6, 7 and 8. Figs. 5 and 7 depicts 
respectively the dependence of the coupling constants <?H e E: c 7r and gs b s b TT on M 2 in the working region of M 2 . The 
results are given for two fixed values of (3 and two fixed values of sq for each coupling constant. It follows from 
the figures that the results are rather stable with respect to the variations of M 2 in the given region of M 2 . The 
dependence of the coupling constants on cos 9 are also presented in Figs. 6 and 8 and with respect to these figures when 
cos# is in between —0.25 < cos 9 < 0.25 the coupling constant 3h c h c it is practically independent of the unphysical 
parameter /3. The interval of /3 that gives coupling constant result independent of f3 for 3H 6 H 6 7r is —0.75 < cos9 < 0.75. 
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As a result of our analysis we obtain the values of the coupling constants as 

SB e s e1 r = 1.0±0.5, «te»s^ = 1.6±0.4. 

To summarize, in this work we present the results of the coupling constant for the coupling constants <7e c e C 7t and 
gs b E b Tr- To this end, we make use of the LCQSR approach with the Sq current applied in its most general form. We 
obtain the appropriate values of the threshold parameters so from the mass sum rules and, using them and appropriate 
intervals of Borel parameter and /3, we attain the coupling constants gs c E c n and <fe b E: b 7r- 
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and /i(of) = V||(ot), ^(a*) = V±(oti), h{cti) = A\\(ai), U{ai) = vA\\(cti), h{ai) = T{ch), /s(ai) = vT(cti), 
hip-i) = •A±(oii) and fs(cti) = vA±(ai) are the pion distribution amplitudes. Note that, in the above equations, the 



are the same, we can set M\ — Mf and u = | 



Borel parameter A/ 2 is defined as A/ 2 = j^s^M 2 an< ^ u ° = m^+m 5 ■ Since the mass of the initial and final baryons 
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FIG. 1: The mass of the H c as a function of the Borel parameter M 2 for different values of arbitrary parameter 
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FIG. 2: The dependence of the mass of the H c on cos 9 for different values of Borel parameter M and the continuum threshold 
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FIG. 3: The mass of the H;, as a function of the Borel parameter M 2 for different values of arbitrary parameter /3 and the 
continuum threshold so- 




FIG. 4: The dependence of the mass of the E b on cos 8 for different values of Borel parameter M 2 and the continuum threshold 
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FIG. 5: The dependence of the coupling constant gs c 
ft and the continuum threshold so- 



on the Borel parameter M 2 for different values arbitrary parameter 
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FIG. 6: The dependence of the coupling constant gs c s C Tr on cost9 for different values of Borel parameter M 2 and the continuum 
threshold so. 
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FIG. 7: The dependence of the coupling constant gs 
P and the continuum threshold so- 



on the Borel parameter M 2 for different values of arbitrary parameter 
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FIG. 8: The dependence of the coupling constant g~. b s. b -n on cos 9 for different values of Borel parameter M and the continuum 
threshold so. 



